Similarity Renormalization Group Evolution of Three-Nucleon Forces 
in a Hyperspherical Momentum Representation 
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A new framework for computing the Similarity Renormalization Group (SRG) evolution of three- 
nucleon forces (3NF) in momentum representation is presented. The use of antisymmetric three- 
particle hyperspherical momentum states ensures unitary evolutions within certain basis truncations, 
much like antisymmetric harmonic oscillator SRG evolutions. Additionally, in each partial wave 
the Tj.gj-SRG regulator is exactly represented, similar to recent 3NF momentum representation 
evolutions. Unitary equivalence is demonstrated for the triton using several chiral two- plus three- 
nucleon interactions. This method allows for a clean visualization of the evolution of the three- 
nucleon forces, which manifests the SRG decoupling pattern and low-momentum universality. 

PACS numbers: 21.30.-x,05.10.Cc,13.75.Cs 



The Similarity Renormalization Group (SRG), as for- 
mulated for nuclei in Refs. [IHSj, renormalizes and 
thereby softens inter-nucleon interactions. One solves a 
flow equation to generate a unitary flow of the Hamilto- 
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where T^q\ is the relative kinetic energy and V consists 
of all inter-nucleon interactions. The details of the flow 
are controlled by choosing an anti-hermitian generator 
rjs ■ Most of the studies in nuclear physics to date use 
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The Tj.gj-SRG flow of two-nucleon forces (2NF) has been 
studied in detail using both momentum representation ^ 
llHl] and in a discrete harmonic oscillator (HO) basis [8]. 
In contrast, three-nucleon forces (3NF) have primarily 
been studied in the HO basis [5HT^ with only recent work 
in momentum representation |13l I14j . 

We have developed an alternative momentum represen- 
tation SRG evolution that exploits hyperspherical har- 
monics (HH) to build a hybrid method combining the 
relative strengths of the previous HO SRG implemen- 
tations and the recent momentum representation SRG 
implementations. This is achieved by representing the 
three-body Hamiltonian in antisymmetric HH momen- 
tum states, which leads to a permutationally closed trun- 
cated representation of the interaction. Our hybrid mo- 
mentum representation method provides evolved interac- 
tion matrix elements that are applicable directly in infi- 
nite systems and are easily projected to a HO basis with 
arbitrary frequency for use in finite systems. Extension 
to alternative SRG generators and four-body interactions 
are straightforward. Finally, many features of the SRG 
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fiow for the 3NF become easy to visualize in a HH mo- 
mentum representation. 

The partial- wave momentum representation fiow equa- 
tion for the two-nucleon system can be written as 



^V^^'^ik, k') = -(fc2 _ fc'')2T/"'^(fc, k') 

as 
+ J2 / I'dq {k' + k'' - 2g2)T/"-7(fc^ q)K'''f^{q, k'), 

(3) 

where Greek letters abbreviate the collection {I, s,j, t, t^) 
of orbital angular momentum, spin, total angular mo- 
mentum, isospin, and isospin projection quantum num- 
bers. The first term on the right side of Eq. (Is]) drives 
the decoupling and also suggests a different parameter to 
use for studying the Tj.gj-SRG: 
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the reduced mass 



for the relative Hamil- 
1 



where /i is 

tonian and A has units of fm~^. When solving the 
T^.g]^-SRG fiow equation for the 2NF, the partial waves 
decouple into 1 x 1 or 2 x 2 blocks, which can be evolved 
separately. This can be done in either Jacobi momentum 
representation with a quadrature rule, or in a discrete ba- 
sis representation such as a truncated set of HO states. 
In either case, the flow is unitary up to the precision of 
the ODE solver used. 

Solving the flow equation for the evolved 3NF using 
an antisymmetric Jacobi HO basis, it is only necessary 
to construct the sum of the initial two- and three-body 
potentials in an antisymmetric basis and evolve using 
Eq. (fTl) with Vs = 3y,^^ + y/^ and the three-body 
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T-^qI (see [5]). Because the antisymmetrizer is block di- 
agonal in this basis, and the initial truncated Hamilto- 
nian is completely contained in some finite set of blocks, 
the SRG fiow equation will not induce non-zero matrix 
elements into blocks outside of the original truncation. 



The evolved 2NF matrix elements can be embedded 
and subtracted from Vg to isolate the evolved three-body 
force, V/^ = Ki-SV/^. While being formally exact, fi- 
nite numerical precision means the subtraction leaves the 
evolved 3NF with some small residual from the evolved 
2NF matrix elements. This may have an effect when us- 
ing the evolved 3NF in a larger nucleus [HKT^]), however 
this effect has not been studied in detail. An additional 
issue with using a discrete basis is that the Tj.g^-SRG flow 
equation induces a regulator that is naturally described 
in a continuum representation. Evolving in momentum 
representation and then embedding in the HO basis can 
yield different evolved matrix elements than evolving di- 
rectly in the HO basis due to truncated contractions in 
the flow equation. Additionally, for each oscillator pa- 
rameter used in an HO basis calculation, the SRG evo- 
lution needs to be recomputed. For these reasons, it is 
more natural to solve for the three-body flow equations 
forces in momentum representation and then project into 
the HO basis. 

In a recent implementation of a momentum represen- 
tation SRG 3NF ^3j , the 3NF is represented in a partial- 
wave basis, 

\kik2;[ihsi)jiil2^)j2]J;{hl)TT,) = \hk2;a), (5) 
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where subscript 1 refers to the relative coordinates for 
particle 1 and 2, subscript 2 refers to relative coordinate 
of the third particle in the center of mass frame of par- 
ticles 1 and 2, and J,T,Tz are total coordinates for the 
three-body systems. Solving the three-body SRG flow 
equation in this momentum representation comes with 
two related difficulties. Expressing the 2NF in this ba- 
sis, there is a spectator delta function that is numeri- 
cally difficult to represent. To overcome this, Ref. [13] 
applied a separated form of the the flow equation from 
Ref. [1 . This form expresses the flow equation as coupled 
flow equations for the 2NF and 3NF, where in the 3NF 
flow equation, the 2NF spectator delta function can inte- 
grated out analytically. Upon solving the flow equation, 
isolated solutions for the 2NF and 3NF are generated 
automatically, and the resulting forces can be projected 
into a HO basis with any oscillator parameter for use 
in many-body calculations. Since the 3NF flow equa- 
tion is explicitly separated from the 2NF flow equation, 
the possible 2NF contamination in the 3NF matrix ele- 
ments, which can occur in the HO SRG evolutions, is not 
present. A second difficulty is that anti-symmetrization 
of the 2NF is not block diagonal within any partial- wave 
truncation of this basis. This enters through terms in 
the flow equation such as AV^^A x Tj.gj x V"^^ , where 
yA is the antisymmetrizer, will induce non-zero matrix 
elements in partial waves outside the initial truncation. 
Ignoring these matrix elements generates a unitary vi- 
olating error. In current evolutions based on |13j . this 
error is at the ^ 0.5 keV level for the triton [IS], but the 



scaling effects of this error in heavier systems are cur- 
rently unknown. 

By exploiting the HH formalism (see [inillT] for details 
of the HH basis), we have developed a hybrid between 
the two approaches for computing the SRG 3NF. We ex- 
press the Hamiltonian using fully antisymmetric Jacobi 
three-body hyperspherical harmonic momentum states. 
Starting from the representation of [13j we project the 
three-body Jacobi momentum states onto hyperspherical 
harmonics. 
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where Q = \/k\ + k\ is the hyper-momentum and G 
is the grand angular momentum. Using results from 
Refs. [T71 [IS], we embed the 2NF between particles 1 
and 2 into the three-body HH momentum basis. 

For identical particles (equal mass), all dependence of 
any orthogonal transformation operator (O), such as per- 
mutation operators, is completely encoded between the 
partial waves and is block diagonal in G [T^, therefore 
they are independent of the remaining continuous coor- 
dinate Q, 

(Q'G';a'|0|QG;a) = (G;a'|0|G;a)JG',G^^^^^. 

(9) 
We build the antisymmetrizer A. in blocks of fixed G, 
diagonalize it, and construct a set of linearly independent 
antisymmetric states in a manner similar to Refs. [20ll21j . 
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We project from our initial basis into a complete set of 
non-degenerate antisymmetric states that have definite 
Q and G, 



{Q'G'i'\V\QGi) 
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This will also build in interactions between particles 1, 3 
and 2, 3. Now the Hamiltonian is translationally invari- 
ant and antisymmetric, computed in the manner of |21j . 
but using a momentum representation instead of a dis- 
crete HO basis. 

This antisymmetric representation can be inserted into 
the SRG flow equation. In this form, the flow equation 
for the total 2-I-3NF is nearly identical in form to Eq. pi, 

/•OO 

+V / P''dP{Q^+Q'^-2P^)V^'^{Q,P)V^'\P,Q'), 

c Jo 

(13) 




FIG. 1: (color online) HH SRG evolution of the triton bind- 
ing energy using a N^LO 500/600 MeV interaction, including 
only the 2NF, the 2NF + induced 3NF, and the 2NF + full 
3NF (see Ref. g]. The dashed lines are for the other N^LO 
2NFs+3NFs( See [^E5|). 



with a, b, c each specifying a hyperspherical channel 
\Gi; JTTz). A critical difference between this flow equa- 
tion and the flow equation solved in |13j is that our 
equation cannot induce non-zero matrix elements above 
our partial-wave truncation, provided all antisymmetric 
states up to some Gmax are included. This enforces 
unitary observables even when using a finite momen- 
tum partial-wave basis. Finally, we isolate the 3NF from 
Vs{Q,Q') in a manner identical to what is done for HO 
evolutions. Full details on our embedding and flow equa- 
tion will be presented in a following work j27j. 

Using the Sundials CVODE solver [21j , we evolved sev- 
eral triton Hamiltonians, setting the proton and neutron 
mass to twice the proton-neutron reduced mass. This is 
done to simplify the evolution and embedding of different 
2NF isospin channels for this initial work. We truncate 
the hyperspherical harmonic basis at Gmax — 20. We 
use the chiral N^LO 2N-h3N Hamiltonians from [22l ES] 
to explore unitary renormalization of the chiral effective 
theory with inital 3NFs while allowing a direct compari- 
son with the 3NF SRG evolutions presented in Ref [T5] . 
These approximations give a result for the triton bind- 
ing energy within ~ 3 keV of the expected energy (fit 
to —8.482 MeV [13]), using 60 momentum states in each 
partial wave. Analyzing our convergence with respect to 
Gmax indicates that the remaining error is dominated by 
our approximate handling of the different nucleon masses, 
and finite quadrature effects. 




FIG. 2: (color online) SRG evolution of the triton probability 
distribution as a function of hypermomentum Q for several 
different SRG As for the interactions in Fig. IT] These are 
plotted as bands that span the range of the wave functions 
from different initial 2N+3N N^LO interactions. 



Figure [T] shows the SRG evolution of the triton binding 
energy for the N'^LO 550/600 MeV interaction using only 
2N and 2N-I-3N induced forces, as well as the binding en- 
ergy for various the N^LO 2NFs + 3NFs. Using only 
the evolved 2N force, we see a non-unitary fiow. This is 
from a missing induced 3NF, which when included yields 
a unitary flow. This behavior is identical to what has 
been documented in Refs. [8} 113). Since we maintain ex- 
act antisymmetry in our finite truncated basis, the error 
in our ground-state eigenvalue is determined by the pre- 
cision of our ODE solver. For the evolutions presented, 
we find this error to be about 1 eV. 

In Fig. [2] we plot the evolution of the triton momen- 
tum distribution, Q^ Tlci {Q^A"^) ■ Each color band is 
the span of the hyper-momentum probability for a set 
of different initial 2N-I-3N N LO interactions. As the 
SRG renormalizes the interactions, we see two critical 
features. The first and most drastic is the suppression 
of the high-momentum tail. This feature has been well 
studied for the dcuteron [3l [S^ and is expected for all 
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■SRG evolved wave functions. The second critical 



feature is that the range spanned by wave functions for 
different interactions collapses to a universal form, signal- 
ing low-momentum universality, again in a manner very 
similar to what has been seen in the deuteron [5S] . This 
collapse is also seen in individual partial waves. 

The interaction is plotted in Figs. [3] and |4] for the low- 
est fully antisymmetric partial wave along the top row 
as a function of the SRG flow parameter A. Along the 
bottom row is a two-body channel, a mixture of ^Sq and 
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FIG. 3: (color online) Color Contour plot of the 3NF (upper row) and 2NF (lower row) as a function 
MeV potential. The lowest antisymmetric H.H. partial wave is plotted as well as the two-body partial 
kernel for for this three-body partial wave. 
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of the A for the 550/600 
waves that in embedded 
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FIG. 4: (color onhne) Same as Figure [s] but using the 600/700 MeV potential. 
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FIG. 5: (color online) Contour plot of the integrand of triton 3NF expectation 

f|*i(Q)*]((5')(W)'^^i?f^(Q>Q')- The upper row is the integrand for several initial N^LO 
the lower row is the integrand for the SRG evolved interactions. 



value, I{Q,Q') = 
chiral potentials, while 



^Si , that is embedded in this three-body partial wave. As 
the interaction is softened via the SRG, we see the band 
diagonahzation/decoupUng that has been described for 
the 2NF in Ref. T], which is as expected from the flow 
equation Eq. ( [l3| . Also simflar to the 2NF SRG results, 
we find that the 3NF interactions are collapsing to a uni- 
versal low-momentum form, though the degree of univer- 
sality is not as extreme as seen for the 2NF evolutions 
(see HI)- We see a similar pattern in the other low 
lying partial waves. 

In Fig. [5j we plot the integrand for triton expectation 
value of the 3NF, 

I{Q:Q')^{Q'Q)' E (*|QGz)(Q'GV|*) 



G,i;G'i' 



x{QGi\V^'^\Q'G'i'). (14) 



The SRG moves critical strength from the high- 
momentum region of the integrand to low momentum. 
In the process, it develops universal integrands for the 
3NF's expectation value, rendering the low-momentum 
universality obvious. 

A feature in our interactions is an apparent non- 
smoothness at small momenta ('^ 1 fm~ ) and near the 
diagonal. These can seen in Fig. |3] near the diagonal 
around 1 fm^ for A = 1.40 fm^ and in Fig.plin the lower 
panels. This error is not seen when projecting matrix el- 
ements from [13) into our representation. Therefore, it 
seems to be a result of computing the 3NF by subtracting 
the evolved 2NF from the evolved 2NF-I-3NF matrix el- 
ements, similar to the HO evolution. This suggests that 
this small error may exist in the oscillator matrix ele- 
ments that have been used in recent calculations (such 
as UH [2S] ) , which could become significant in larger 
systems and should be studied further. 

Another advantage of our approach is that our mo- 
mentum representation Hamiltonian can be completely 



diagonalized on a given momentum quadrature. Exploit- 
ing this, is it possible to construct the unitary operator 
describes the SRG, permitting the evolution of various 
operators within momentum space. 

We have developed a fully unitary momentum repre- 
sentation SRG 3NF evolution and applied it to chiral 
N LO inter-nucleon interactions in the triton channel. 
This was accomplished by expressing the Hamiltonian in 
fully antisymmetric hyperspherical plane waves. Exam- 
ining the evolution in this representation, we find that all 
features observed in 2NF momentum representation evo- 
lutions are manifest in the triton system: suppression of 
high-momentum tail, collapse of wave functions and in- 
teractions from different initial Hamiltonians to univer- 
sal forms, decoupling between high- and low-momentum 
states, unitarity of the eigenvalues to the precision of 
ODE solver used. In the process, we identified features 
of evolved 3NF interactions that are a result of the finite 
precision cancellation of the 2NF, for which the nature 
and effect in heavier systems needs to be investigated. 
Using the HH momentum representation, the momentum 
representation evolution of the four-nucleon force should 
be straightforward. 

We gratefully acknowledge R. J. Furnstahl, K. Hebeler, 
and H. Hergert for many helpful discussions, and R. 
J. Perry, S. More and B. S. Dainton for useful com- 
ments, and E. Epelbaum and A. Nogga for provid- 
ing the initial N LO 3NF matrix elements to the OSU 
group and N. Barnea for notes on HH calculations. 
I thank the Ohio Supercomputing Center for the use 
of their computational resources. This work was sup- 
ported in part by the National Science Foundation un- 
der Grants No. PHY-0758125 and No. PHY-1002478, 
the UNEDF SciDAC Collaboration under DOE Grant 
DE-FC02-07ER41457, and by the DOE Office of Sci- 
ence Graduate Fellowship (SCGF) program under DOE 
contract number DE-AC05-06OR23100. 



S. K. Bogner, R. J. Furnstahl, and R. J. Perry, Phys. 

Rev. C 75 (2007). 

S. K. Bogner, R. J. Furnstahl, R. J. Perry, and 

A. Schwenk, Phys. Lett. B 649, 488 (2007). 

E. D. Jurgenson, S. K. Bogner, R. J. Furnstahl, and R. J. 

Perry, Phys. Rev. C 78, 014003 (2008). 

K. A. Wendt, R. J. Furnstahl, and R. J. Perry, Phys. 

Rev. C 83 (2011). 

K. A. Wendt, R. J. Furnstahl, and S. Ramanan, Phys. 

Rev. C 86, 014003 (2012). 

W. Li, E. R. Anderson, and R. J. Furnstahl, Phys. Rev. 

C 84, 054002 (2011). 

E. R. Anderson, S. K. Bogner, R. J. Furnstahl, E. D. 

Jurgenson, R. J. Perry, and A. Schwenk, Phys. Rev. C 

77, 037001 (2008). 

E. D. Jurgenson, P. Navratil, and R. J. Furnstahl, Phys. 

Rev. Lett. 103, 082501 (2009). 

E. Jurgenson, P. Maris, R. Furnstahl, P. Navratil, W. Or- 



mand, et al. (2013), 1302.5473. 
[10] E. D. Jurgenson, P. Navratil, and R. J. Furnstahl, Phys. 

Rev. C 83, 034301 (2011). 
[11] S. Binder, J. Langhammer, A. Calci, P. Navratil, and 

R. Roth, Phys. Rev. C 87, 021303 (2013), ISSN 0556- 

2813. 
[12] R. Roth, S. Binder, K. Vobig, A. Calci, J. Langhammer, 

and P. Navratil, Phys. Rev. Lett. 109, 052501 (2012). 
[13] K. Hebeler, Phys. Rev. C 85, 021002 (2012). 
[14] O. Akerlund, E. J. Lindgren, J. Bergsten, B. Grevholm, 

P. Lerner, R. Linscott, C. Forssen, and L. Platter, Eur. 

Phys. J. A 47, 122 (2011). 
[15] K. Hebeler, private communicaton. 
[16] N. Barnea, W. Leidemann, and G. Orlandini, Nucl. Phys. 

A 693, 565 (2001). 
[17] M. Viviani, L. E. Marcucci, S. Rosati, A. Kievsky, and 

L. Girlanda, Few-Body Syst. 39, 159 (2006). 
[18] A. Novoselsky and N. Barnea, Phys. Rev. A 51, 2777 



(1995). 
[19] J. Raynal and J. Revai, Nuovo Cimento A 68, 612 (1970). 
[20] N. Barnea, Phys. Rev. A 59, 1135 (1999). 
[21] P. Navratil, G. Kamuntavicius, and B. Barrett, Phys. 

Rev. C 61 (2000). 
[22] E. Epelbaum, Progress in Particle and Nuclear Physics 

57, 654 (2006), ISSN 01466410. 
[23] E. Epelbaum, A. Nogga, W. Glockle, H. Kamada, U.-G. 

Meifi ner, and H. Witaa, Phys. Rev. C 66, 064001 (2002). 
[24] A. C. Hindmarsh, P. N. Brown, K. E. Grant, S. L. Lee, 



R. Serban, D. E. Shumaker, and C. S. Woodward, ACM 
Transactions on Mathematical Software 31, 363 (2005). 

[25] E. R. Anderson, S. K. Bogner, R. J. Furnstahl, and R. J. 
Perry, Phys. Rev. C 82, 054001 (2010). 

[26] H. Hergert, S. K. Bogner, S. Binder, a. Calci, J. Lang- 
hammer, R. Roth, and a. Schwenk, Phys. Rev. C 87, 
034307 (2013). 

[27] K. A. Wendt, R. J. Furnstahl, in preparation. 



